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The estimation of cosmological parameters from cosmic microwave experiments has almost always 
been performed assuming gaussian data. In this paper the sensitivity of the parameter estimation 
to different assumptions on the probability distribution of the fluctuations is tested. Specifically, 
adopting the Edgeworth expansion, I show how the cosmological parameters depend on the skewness 
of the Ci spectrum. In the particular case of skewness independent of £ I find that the primordial 
slope, the baryon density and the cosmological constant increase with the skewness. 
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I. INTRODUCTION 

The new generation of cosmic microwave background 
(CMB) experiments (see e.g. fil, g, |J, and the future 
missions Planck [Q, and Map flBj) promises to estimate 
the cosmological parameters within a precision of 1%. 
The current dataset already allows in some cases an un- 
certainty below 10% on such parameters as the baryon 
density or the primordial spectral slope. Such a precision 
allows and demands a clear assessment of the theoretical 
assumptions. 

So far, all the estimations of cosmological parameters 
based on the CMB data assumed a Gaussian distribu- 
tion of the primordial fluctuations (the only exception I 
know of is Ref. || in which the primordial slope n was 
estimated in presence of skewness). Such an assump- 
tion is based on the conventional models of inflation and 
has the obvious and enormous advantage of being sim- 
ple and uniquely determined. However, the gaussianity of 
the primordial temperature fluctuations is still to be fully 
tested |0, ||, y and there exist several theoretical models 
which actually predict its violation |1C|| . Therefore, it is 
necessary to quantify how the cosmological parameters 
derived from CMB experiments depend on the statistical 
properties of the fluctuation field. 

In this paper I derive the dependence of four cosmo- 
logical parameters on the skewness of the fluctuations as- 
suming a flat space. The parameters are the primordial 
slope n, the baryon and cold dark matter rescaled density 
parameters u>b = Q.bh 2 ,uj c = Q c h 2 ( h is the Hubble con- 
stant in units of 100 km/sec/Mpc), and the cosmological 
constant density parameter SI a. The flat space constraint 
reduces to the relation h 2 = (w& + W c )/(1 — 0,\). 

It is clear that removing the hypothesis of Gaussianity 
leaves room for an infinity of different possible assump- 
tions concerning the fluctuation distribution. I adopt 
here the Edgeworth expansion (EE), for three reasons: 
a) it can be seen as a perturbation of a Gaussian func- 
tion; b) it is easy to manipulate analytically and c) it is 
the distribution followed by any random variable that 
is a linear combination of N random variables in the 
limit of large N (for N — » 00 the Edgeworth distribu- 
tion reduces to a Gaussian). The latter property might 
be useful to describe fluctuations that arise due to several 



independent sources. The EE has been previously used 
in cosmology to model small deviations from Gaussianity 

The main drawback of the EE is that it is not positive 
definite. However, when the deviation from Gaussianity 
is small, this problem is pushed many standard deviations 
away from the peak and does not affect the parameter 
estimation. 

This paper is meant to exemplify the effects that a non- 
zero skewness introduces on the likelihood estimation. 
For generality, I will not confine myself to any specific 
mechanism for generating the non-gaussianity. Moreover, 
for simplicity, I will skip over several additional compli- 
cations like bin cross-correlations, calibration, pointing 
and beam errors that an accurate analysis should take 
into account. 



II, EDGEWORTH LIKELIHOOD 

The likelihood function usually adopted in CMB stud- 
ies (e.g. [Q, g, |^, [16|) is an offset log-normal function. 
This function is an approximation to the exact likelihood 
that holds for Gaussian data in presence of Gaussian 
noise | fi6| . The offset depends on quantities that are not 
yet publicly available; since the offset can be neglected 
in the limit of small noise, we assume as starting point 
a simple log-normal that, neglecting factors independent 
of the variables, can be written as 



21ogL( aj )=J2 



Zt,t(2ii a j) ~ Z l,d(£i) 



(1) 



where Z? = log Cg, the subscripts t and d refer to the the- 
oretical quantity and to the real data, Cg are the spectra 
binned over some interval of multipoles centered on £i , at 
are the experimental errors on Zt t d> and the parameters 
are denoted collectively as a? . We neglect also the resid- 
ual correlation between multipole bins, which should be 
anyway very small for the latest data. An overall am- 
plitude parameter A can be integrated out analytically 
adopting a logarithmic measure dlogA in the likelihood. 
Writing Ci = AC' e it follows Z iyt = log A + log C' t t = 



log A + Z' g t so that, neglecting the factors independent 
of the variables and putting w = log A, we obtain 
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Let us now introduce the Edgeworth expansion. De- 
noting with xt = (w + At) /at the normal variable in the 
Gaussian function, the Edgeworth expansion is p^ 



P(xi) = exp 
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where k n> i is the n-th cumulant of xt and H n is the Her- 
mite polynomial of n-th order. Notice that the EE has 
the same norm, mean and variance as the Gaussian, but 
different mode (the peak of the distribution). Here, as 
a first step, we limit ourselves to the first non-Gaussian 
term containing the skewness ks,t- 

Assuming that X£ is distributed according to the Edge- 
worth expansion, we can build the truncated Edgeworth 
likelihood function pi] to first order in fc 3 ^: 



= « 2se 
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with H%{xi) — x\, — ixt- Now, integrating over w we 
obtain 
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where 
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This is the likelihood function that we study below. The 
effect of the extra terms is to shift the peak (or mode) of 
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Figure 1: The squares mark the mode of the G( distribution 
for &3 = —2 (light squares), and for ks = 2 (dark squares). 
The data are from Cobe and Boomerang. 



the distribution of each Ct while leaving the mean unper- 
turbed. Since the shift depends on At ,at and k3,t, the 
resulting mode spectrum will be distorted with respect 
to the mean spectrum. Therefore, the likelihood maxi- 
mization will produce in general results that depend on 
k 3 ^- In Fig. 1 we show the peak shift introduced in the 
simplified case in which the skewness is independent of 
£: if k 3 is negative, the spectrum is shifted upward by 
a larger amount at the very small and very large multi- 
poles, and by a smaller amount around £ — 200, where 
the relative errors are the smallest; if k 3 is positive the 
shift is downward. As a consequence of the distortion, we 
expect that a constant negative skewness favours spectra 
which are tilted downward with respect to the Gaussian 
case, and the contrary for a positive skewness. In general, 
the cosmological parameters will depend on the multipole 
dependence of k 3< £. For small k 3 at, the shift can be ap- 
proximated by 



Ct(mode) — CV( mea „)(l — k 3 at/2). 



(5) 



Clearly, if the peak shift introduced by the EE were inde- 
pendent of £, the integration over the amplitude w would 
erase the non-gaussian effect on the likelihood. That is, 
putting at and A^ equal to a constant independent of £ 
we obtain g(k 3t i, At, at) = 0. 



III. DEPENDENCE ON THE SKEWNESS 

To evaluate the likelihood, a lib rary of CMB spec- 
tra is generated using CMBFAST @. Following j| 
I adopt the following uniform priors: n E (0.7,1.3), 
Lu b E (0.0025,0.08), uj c E (0.05, 0.4), fl A E (0,0.9). 
As extra priors, the value of h is confined in the range 
(0.45,0.9) and the universe age is limited to > 10 Gyr. 
The remaining input parameters requested by the CMB- 
FAST code are set as follows: T cmb — 2.72QK, Y He = 
0.24, N v = 3.04, t c = 0. In the analysis of r c , the op- 
tical depth to Thomson scattering, was also included in 



the general likelihood and, in the flat case, was found 
to be compatible with zero at slightly more than la . 
Therefore here, to reduce the parameter space, I assume 
t c to vanish. The theoretical spectra are compared to the 
data from COBE fil| and Boomerang @. 

To specify the skewness k 3t g three simplified cases are 
studied: in the first one ("constant skewness"), k 3 ,g = k' 3 
is assumed independent of the multipole £; in the sec- 
ond ("gaussian skewness") , the skewness is assumed to be 
generated by some process only in a particular range of 
multipoles: 
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where, in the numerical examples below, I put £* = 200 
and Ag = 20. In the third case ("hierarchical skewness"), 
the "hierarchical" ansatz is assumed Q18[ [ , in which the 
skewness of the temperature field is proportional to the 
square of its variance. At the first order, we can assume 
that the skewness of the Cg distribution is proportional 
to the skewness of the fluctuation field, so I put 



k 3 , e = k; {Ct/cry 



(7) 



where, for instance, C* — 6*200- In all three cases fcj is 
left as a free parameter. These three choices are of course 
purely an illustration of what a real physical mechanism 
might possibly produce. 

Fig. 2 shows the one-dimensional Edgeworth likeli- 
hood functions marginalized in turn over the other three 
parameters. For the "constant skewness", k 3 varies from 
-1.6 to 1.2 ( light to dark curves): below and above these 
values the likelihood begins to show pronounced nega- 
tive wings, which signals that the first order Edgeworth 
expansion is no longer acceptable. While the likelihood 
for lo c is almost independent of k 3 , it turns out that the 
other likelihoods move toward higher values for higher 
skewness. As anticipated, this can be explained by ob- 
serving that a higher skewness implies smaller Cg at small 
multipoles: a tilt toward higher n and higher u>b gives 
therefore a better fit. The effect is of the order of 10% 
for fc^ w 1. 

In the "gaussian skewness" case the trend is qualita- 
tively the opposite, as can be seen in Fig. 3, where /c| 
ranges from -4 to 4 (light to dark). Here the cosmo- 
logical parameters decrease for an increasing skewness. 
The reason is that now the effect is concentrated around 
the intermediate multipoles £ ~ 200: a positive skewness 
induces smaller Cg at these multipoles, and therefore a 
smaller n and ujb helps the fit. The third case, the "hier- 
archical skewness", is not shown because is qualitatively 
similar to the previous case: the region around £ = 200 
is in fact also the region where Cg is larger and therefore 
k3 t t given by Eq. (0) is larger. 

Fig. 4 summarizes the results: the trend of the esti- 
mated parameters (mean and standard deviation) versus 
feg in the "constant skewness" case. The constant plateau 
that is reached for k% < depends on the fact that for 
large and negative skewness the peak shift is independent 
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Figure 2: Likelihood functions for the four cosmological pa- 
rameters in the constant skewness case. The skewness in- 
creases from -1.6 to 1.2, light to dark. 




Figure 3: Likelihood functions for the four cosmological pa- 
rameters in the "gaussian skewness" case. The skewness in- 
creases from -4 to 4, light to dark. 



of &3 . The cosmological parameters can be well fitted by 
the following expressions: 



n= 0.90(1 + 0.03e fe3 ), 

uj b = 0.021(1 + 0.05e fe3 ), 

w c = 0.115(1 + 0.025e fe3 ), 

Qa= 0.67(1 + 0.06e fc3 ). 



(«) 

(9) 

(10) 

(11) 



For h the fit is h 2 = 0.42 ^+°-°g e £ . Notice that the trend 

1— 0.12e"-3 

for h is stronger than for the other variables: h goes from 
0.65 to 0.85 when & 3 increases from -1.6 to 1.2. Similar 
relations can be found for the other cases as well. 




Figure 4: Variation of the mean and standard deviation of the 
cosmological parameters versus k%, in the "constant skewness" 
case. 



IV. CONCLUSIONS 

This paper illustrates quantitatively a basic and ob- 
vious fact about cosmological parameter estimation, 
namely the dependence on the underlying statistics. Al- 
though the gaussianity of the CMB data is still to be 
proved, almost all the previous works estimated the cos- 
mological parameters assuming vanishing higher order 
cumulants. Here it has been shown that a non-zero skew- 



ness distorts the mode spectrum with respect to the mean 
spectrum, inducing a considerable variation to the best 
fit cosmological parameters. 

The Edgeworth expansion we used in this paper is con- 
venient for analytical purposes but its use is limited to 
relatively small deviations from gaussianity. In fact, the 
peak shift displayed in Fig.l is always smaller than the 
errobars, and as a result the parameters, although vary- 
ing with ka, remain always within one sigma from the 
zero-skewness case. This, however, does not mean that 
the dependence on the higher order moments can be ne- 
glected, first because it is a systematic effect, and second 
because more general probability distributions which are 
not small deviations from gaussianity might introduce 
much larger shifts. 

We have shown that, to first order, the peak shift 
AC'i/Ci is proportional to k^^ai. The error at includes 
cosmic variance and experimental errors. In the future, 
the main source of error will be cosmic variance, at least 
below I — 2000 or so. A skewness of order unity will 
therefore introduce an additional "skewness bias" that 
will limit the knowledge of the cosmological parameters 
by an amount similar to the cosmic variance itself. At 
this point it will become necessary to estimate k^t along 
with the other parameters. The first order EE is however 
inadequate, since it is linear in k 3 j, and it will be neces- 
sary to extend the expansion to higher orders || [U|, or 
to adopt a non-perturbative non-gaussian distribution. 
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